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Abstract 

We study the prospects for observing double parton scattering through four- 
muon final states, forming two opposite-sign muon pairs, in the LHCb experi- 
ment in pp collisions at 14 TeV centre of mass energy. We consider two special 
cases, namely double Drell-Yan and J/^/^-pair production. The kinematic prop- 
erties and prospects for observing these processes are discussed. We find that the 
production rate depends strongly on the origin of the four muons, while many 
kinematic properties can be used to help identify the presence of double parton 
scattering events. 

1 Introduction 

The subject of multiple parton interactions has been enjoying a renewed interest over 
the last few years, driven by the need to understand the full spectrum of hadronic 
activity at the LHC The large collision energy available at the LHC increases the 
probability of scattering of more than just one parton from each hadron. In particular, 
multiple scattering events with final state particles carrying relatively high transverse 
momentum should occur much more commonly than at the previously operating hadron 
colliders. Such observations provide valuable information on the structure of the proton 
and the parton-parton correlations within. Indeed, processes involving two hard scat- 
terings in one hadron-hadron collision, the so-called double parton scattering (DPS), 



have been studied in the 4-jet and the 7+3-jet channels by the AFS [Ij, UA2 [2[, 
CDF jif and DO j3] experiments, and the presence of DPS is now well established. 

In order to measure double parton scattering, one needs to consider processes for 
which the single parton scattering (SPS) background is low or possesses characteristics 
that allow for the signal to be extracted. To date, phenomenological studies have 
focused on observation of DPS in hadronic processes such as fourjets |5l-ll3| , bb pairs 14 , 



photons [18|, [19 




15| . bb pair and two jet s |lOl| . Higgs and a gauge boson [16|, ll7|, production of either 



or Z bosons [22[ in association with jets, or purely lepton 



events such as double Drell-Yan [23|, UM same-sign W-pair production [15|, |22|, |25 



26 



In the environment of hadronic collisions, purely leptonic signatures offer a clean 
probe of the underlying scattering mechanisms. A generic property of DPS processes 
is that they are expected to peak strongly in the low transverse momentum (pr), 
low phase space region, which is difficult to access with a typical central detector. 
In final states involving jets, tagging and identifying low px jets associated with the 
DPS signal process is challenging due to the overwhelming QCD background. However 
multi-lepton, especially muon, final states are much cleaner, making triggering on multi 
lepton final states with low pt possible. In this paper we take advantage of the excellent 
low Pt muon acceptance of the LHCb detector, which can go down to ~ 1 GeV, and 
focus on four-muon final states that form two opposite-sign (OS) muon pairs in this 
particular experiment. This should maximize the number of DPS events observed. 

Depending on the analysis cuts, in particular on the invariant mass of the OS muon 
pairs, the DPS signal could be dominated by double Drell-Yan (DDY) production or 
production of a quarkonium (e.g. J/ip or T) pair, which then decay into four leptons. 
The possibility of observing DPS in pair-production of quarkonia has been recently 
investigated in 



27H29I]. In particular, in Ref. [27[ we have argued that the recent 



measurement of double J/tp production by the LHCb experiment S^l could already 
indicate a significant contribution from DPS. 

While the production of double J/ip pairs benefits from a large cross section, the 
subprocess of DDY, i.e. single DY production, is theoretically well understood. For 
this reason, double Drell-Yan production might be considered a standard-candle DPS 
process and merits a more detailed investigation. Moreover, the partons involved in 
these processes are different. DPS double J/ip production is expected to be dominantly 
produced by four gluons, whereas for DDY, the initial states at leading order (LO) are 
two quark-anti-quark pairs. Hence the correlations being probed are different. In this 
paper, we shall primarily focus on the prospect for observing DPS through the DDY 
channel. We also comment on some characteristic differences between double J/ ip and 
DDY production. 

The paper is structured as follows. In Section [2] we discuss the DPS signal and SPS 
background to DDY production at the LHC. Our method of performing calculations 
and cuts chosen to select the signal appropriate to the LHCb experiment are presented 
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Figure 1: Example Feynman diagram for the DPS double Drell-Yan process. 



in Section [3l where we also give the total rates for the processes of interest. In Section 
Hlwe address the differential distributions for the signal and background to the DDY 
process, followed by a discussion of differences between the DDY and double J/ip 
distributions in Section [5l before concluding in Section El 



2 Double Drell-Yan at LHCb 



2.1 DPS signal 

The DPS signal process leading to the two OS muon pair final state is assumed to come 
from two independent hard scatterings of the DY type. An example Feynman diagram 
is shown in Fig. [H If we further assume that the longitudinal and the transverse 
components of the generalized double parton distribution functions can be factorised, 
and that there are no longitudinal momentum correlations between the partons in the 
same hadron, the DPS DDY cross section factorizes into a product of two SPS single 
DY cross sections. More precisely, the differential DPS cross section da^^ is written 
as 



da 



J^DY ^ DY 
DDY "° SPS "'^ SPS 



DPS 



where 



^"^SPS = X] I^F)fb{Xb, IJ'F) da^ps dXadXb : 



(2) 



a,b 
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and ^ 

d^sps = X] 2? \^^b-,r~*i+i-+x\'^ dPSi+i-+x (3) 

a,b 

correspond to the hadronic and partonic single DY process respectively. In the above 
expressions, s is the partonic centre-of-mass energy squared, fi and Xi {i = a, h) are the 
parton distribution functions (PDFs) and longitudinal momentum fraction for parton 
with flavour i, and /i^ is the factorization scale. The single DY matrix element and 
phase space element are denoted by M and dPS respectively. The factor (Jcfr captures 
all the information on the transverse structure of the proton. It is likely to be process 
and energy dependent, and is the main quantity that a DPS experiment will aim to 
extract. However, for concreteness, we shall use a^s = 14.5 mb, a value measured in 
the 7+3-jet study by the CDF experiment |3 |Pj 

Recently, questions have been raised 3l|, |32| about to what extent this intuitive pic- 
ture can be derived in perturbative QCD. In particular, the assumption of factorization 
between the longitudinal and transverse components of double parton distributions ap- 
pears not to be strictly consistent with QCD. It leads to an appearance of the spurious 
"single-feed" term in the evolution equation for double parton distributions, describing 
correlations in longitudinal momentum fractions Xj. However, in the current analysis 
we consider final states produced through DPS with relatively low invariant masses, 
meaning the incoming partons have on average low longitudinal momentum. Since at 
low X the parton-parton correlations are expected to be negligible, using single parton 
distributions fi{xi,fiF) is likely to be a good approximation. 



2.2 SPS background 

Example Feynman diagrams for the SPS production process of two OS muon pairs 
are shown in Fig. [2J They can be classified into two groups: the "double resonance" - 
type diagrams (Fig. [2^), where both virtual photons are attached to the initial state 
quarks, and the "single resonance" -type diagrams (Fig. |2]d), where the quarks fuse into 
a 7* before decaying into four muons via an additional 7*. In our analysis we consider 
contributions from both types of diagrams, as opposed to a previous study |2J] in which 
only the first type of diagrams is included^ 



2.3 Selection of the DPS signal 

As mentioned in the Introduction, in order to observe the DPS signal, the SPS back- 
ground has to be low, at least in some region of the phase space. In the case of two OS 
muon pair final states, there is no particular reason why a priori the SPS production 

^ A similar value, 15.1 mb, was obtained by the corresponding DO experiment^. 
^From (electromagnetic) charge consideration, it can be seen that these two types of diagrams are 
separately gauge invariant. In this sense they could be considered independently. 
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(a) (b) 



Figure 2: Example Feynman diagrams for the SPS double Drell-Yan process, (a) 
"double resonance" -type diagrams, where both virtual photons are attached to the 
initial state quarks, and (b) "single resonance" -type diagrams, where only one virtual 
photon is attached to the quark line. 



rates should be suppressed, as opposed to, e.g., the same-sign TV-pair production, 
where due to electromagnetic charge conservation the SPS process can only occur at 
higher order in powers of (gauge) couplings than in DPS. Therefore in studying the 
final state with two OS muon pairs, it is important to exploit different kinematical 
properties of the SPS and DPS final states. 

For the DPS signal, given our independent hard scattering assumption, each of 
the two OS muon pairs originating from the two 7* will exhibit pt balancing at the 
parton level. The presence of initial state radiation (ISR) will provide a non-zero 
Pt to each of the 7*. The independent scattering hypothesis then implies that the 
azimuthal angular separation between the two 7*s should be uniformly distributed. On 
the contrary, for the SPS background, at the parton level the transverse momentum 
is generally only balanced among all four muons. An implication is that for the two 
combinatoric ways to form two OS muon pairs, which we shall also call 7*, both will 
lead to a configuration in which the two 7*s travel back-to-back in the transverse 
plane. However again this picture will be distorted in the presence of ISR. These 
basic kinematic differences form the basis for identifying the presence of DPS events in 
previous experimental studies. A number of "pair wise balancing" variables, designed 
to exploit the different transverse plane behaviour for the DPS and SPS processes, 
have been introduced in the literature. We shall see that the impact of significant 
QCD radiation on the low invariant masses being considered will lead to significant 
deviation from the intuitive parton-level picture, making the usage of such pair wise 
balancing variables less effective in discriminating DPS from SPS events. On the 
other hand, "longitudinal" variables are generally less sensitive to these effects. The 
absolute rapidity difference between two reconstructed J / ips was first proposed in 27 



as a method to extract DPS from SPS events. We shall see whether a similar strategy 
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can be used for DDY by constructing a similar variable for the 7*s. 

Another complication is that for DDY there is a two-fold (combinatoric) ambiguity 
in grouping the four muons into two OS pairs. This ambiguity is absent for double 
J I'lp, since the "correct" pairing must have fi~^fi~ invariant masses lying within a mass 
window close to the physical J/ ip mass. This leads to a different set of cuts that must 
be applied to the DDY, which will result in different kinematical distributions. We 
shall discuss the impact of these effects in more detail in the following sections. 

So far, we have considered the irreducible SPS background only. In the low mass 
region, semi-leptonic decay of heavy quarks and hadron mis-identification can fake 
the signal, and are important sources of background. These backgrounds have been 
studied in LHCb in the context of single DY 33|] , and it was found that a multi-variate 
selection strategy based on the isolated nature of the muons produced in DY processes 
can be used to achieve purities of 95% down to m^+^~ = 2.5 GeV. A similar method 
could presumably be used for DDY, however a detailed study is beyond the scope of 
the present study. 



3 Event simulation 

We now discuss how the DPS signal and SPS background events are generated. The 
predictions for the DPS signal are obtained using the multi-parton scattering model 



in Herwig++ v2.4.2 [3J|, which allows for generation of two hard scatterings in every 



event. We use the following (default Herwig++ ) a(Q^) definition: a{Q'^) = a(0)/(l — 
K{Q^)), where a(0) = 1/137.04, and K{Q^) = ^(13.4955 + 31og(gVGeV2)) + 
0.00165 + 0.002991og(1.0 + Q^/GeV^). We use MSTW^'oS NLO PDFs |35|, and set the 
factorization scale jj^p for both 2-to-2 single DY subprocesses at jj^p = V^, which can 
be different for the two subprocesses in the same event. Herwig++ enables us to take 
higher order QCD effects into account in the parton shower. In principle Herwig++ can 
also include effects of intrinsic p^^smearing of incoming partons. However currently 
p^-smearing is applied only to the "primary" hard process. In order to provide intrinsic 
p^-smearings for both hard subprocesses, we switch off the default intrinsic Gaussian 
Pt functionality during event generation, and implement independent pt smearings 
to the two hard subprocesses using the same Gaussian model. For illustration, the 
root-mean-square intrinsic pt in this model, a, is set to a conservative value of 2 GeV. 
The SPS events with two OS muon pairs are generated using the package MADGRAPH 



v5. 1.2.4 361]. Here the same a(Q^) definition as in the DPS process is used. We 
include both "double resonance" and "single resonance" -type diagrams, as discussed 
in Section [2l In the low invariant mass region being considered in this analysis, it 
is sufficient to consider only the Feynman diagrams with two 7*, i.e. no diagrams 
with (virtual) Z-bosons are included. For consistency, in the experimental analysis, 
the contributions from processes involving Z's have to be eliminated. This can be 
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achieved by imposing an upper cut on the invariant mass of the four muons, m4^. 
Such a cut is also useful in enhancing the signal-to-background (S/B) ratio, as the 
DPS signal is expected to have comparatively low four-muon invariant mass. The 
factorization scale is set at fip = i.e. 777,4^. The parton level events generated by 
MADGRAPH are interfaced to Herwig++ for parton showering. As in the DPS process, the 
MSTW 08 NLO PDFs are used and an intrinsic Gaussian px smearing is applied after 
event generation, with the parameter a set to 2 GeV. 

We now discuss event selection. For the muons to be observed in the LHCb detector, 
the pseudorapidity (77) and pt of each muon has to obey: 

• 1.9 <r] < 4.9, 

• Pt > I GeV . 

Apart from the muons coming from the decay of 7*, the simulated events can also con- 
tain additional muons generated by the decay of hadrons formed in the hadronisation 
process. To be conservative, we select events with exactly four muons passing the above 
acceptance cuts, and that they form two OS muon pairs. Additionally, we request that 
the events correspond to processes with 7* decaying into muon pairs by avoiding the 
low mass hadronic resonance and the Z resonance regions with the following cuts: 

• i^fi+fj.- > 4 GeV for all four fi'^fi" combinations , 

• < 40 GeV . 

Note that the cut m^+^- > 4 GeV must apply to all jX^ combinations, instead of only 
two of them. This is because there are potentially dangerous backgrounds involving 
low mass hadronic resonances, for example J /ip. These resonances could lead to high 
invariant masses from mismatched pairs, i.e. pairs that do not come from the 

decay of the same resonance. A cut on all fi^fi" pairs will however also apply to the 
true pairing, which of course has a low fi^fi~ invariant mass, and reject the event as a 
result. 

Furthermore, contributions from T(15') and its higher resonances T{2S) and T(35'), 
when combined with other source of muons, can lead to two OS muon pairs. The T 
background can be suppressed by rejecting events with 

• 9.2 GeV < m^+^- < 10.5 GeV for any of the four fi~^fi~ combinations, 

where the mass window is obtained from the recently published result by LHCb on T 



production [37 



We will refer to the cuts above as the set of "basic cuts" for the DDY analysis. 
We assume 100% detection efficiency for an event passing the above cuts. The cross 
sections for the production of four-muon final states in pp collisions at 7 and 14 TeV, 
after applying basic cuts, are displayed in Table [TJ 
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DDY cross sections [fb] at LHCb 




DPS 


SPS 


7 TeV 
14 TeV 


0.08 
0.16 


0.43 
0.68 



Table 1: DPS and SPS DDY cross sections (in fb) for pp collisions at 7 and 14 TeV. 
For clarity, only the results with parton shower and intrinsic Gaussian px broadening 
a = 2 GeV are shown. 



double J/ip cross sections [pb] at LHCb 




DPS 


SPS 


7 TeV 
14 TeV 


3.16 
7.69 


1.70 
2.62 



Table 2: DPS and SPS double J/ip cross sections (in pb), including the branching ratio 
BTl{J/ip — >■ /i"'"/i~). For clarity, only the results with parton shower and intrinsic px 
broadening a = 2 GeV are shown. 



For comparison, we also present predictions for four-muon final states originating 
from a decay of a J/ ?/;-pair, produced in pp collisions at 7 and 14 TeV. Our line of 



analysis follows the one of (27|, where we studied the double J/ip production at LHCb 



for a collision energy of 7 TeV, i.e. we apply the same methods of calculations and the 



same parameter values as in [27[. As for the DDY analysis, the set of single muon cut 
is given by: 

• 1.9 <r] < 4.9, 

• Pt > I GeV . 



The narrower range 2 < 77 < 4.5 in the analysis of |27[ was motivated by the cuts use d by 
the LHCb collaboration in the first experimental analysis of double J/ tp production 30 . 
Out of the two combinatoric ways to form two OS pairs, the combination with the two 
invariant masses closest to the physical J ftp masses is chosen. This requirement is 
sufficient to discriminate between the combinations, as the OS muon pairs originate 
from J/ijj resonances. In our event simulation, we do not include effects due to finite 
detector resolution. As for the DDY analysis, MSTW 08 NLO PDFs are used and 100% 
reconstruction and detection efficiency is assumed. The DPS and SPS cross sections 
obtained after basic cuts and the invariant-mass selection criterion are applied are 
displayed in Table HI 

At present, it is not yet understood how to systematically include higher order ef- 
fects within the DPS framework. The inclusion of the parton shower in our simulations 
may be considered as a possible first step towards predictions with higher order cor- 
rections taken into account. However, the shower evolution does not change the total 
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event rate of the parton level, if integrated over the full phase space. Therefore it is 
illustrative to consider the size of the next-to-leading (NLO) i^-factors for the total 
cross sections. 

In the context of our calculations, the appropriate K-factor would be defined by 
the ratio of the total cross section at NLO calculated using NLO pdfs to the LO total 
cross section using NLO pdfs. For the pp 7*7* — e~^e~fi~^^~ process at 14 TeV, 
calculations with the MCFM package (38| give, for the set of basic cuts considered here, 
Knlo = 1.7. 

For comparison, the K-factor for single DY, defined in the same way as above, is 
-^NLO = 3.0, when the cross sections are calculated with invariant mass m^+^~ between 
4 and 40 GeV. This could indicate an enhanced S/B ratio if higher order effects were 
taken into account. 

Based on the DPS and SPS rates quoted in Tables [T]and [2l we see that the measure- 
ment of DPS in four-muon final states will be more challenging for the DDY process, 
when compared to double J/ijj production. In particular, the rates for DPS DDY are 
of 0(0. 1) fb and 0(0.2) fb at the 7 and 14 TeV LHC respectively, with relatively low 
S/B ratios. However possible improvements could be made. For example, it might be 
possible to better understand the T background by investigating the behaviour of the 
inclusive T sample 37| that contains an extra pair of muons when the "basic cuts" are 
applied. 

We also recall that in our calculation we assume a fixed value for = 14.5 mb. 
This value might be process and scale dependent, and as discussed in Section [2] is a 
parameter to be determined by experiment. Therefore in the next section we will focus 
on studying kinematical distributions and investigate if further kinematical cuts can 
be introduced to improve the S/B ratio, as well as facilitating the extraction of acs- 
With luminosity of a few inverse femtobarn, the number of DDY events expected at 
the 7 TeV LHC is likely to be too low. For this reason only results for the 14 TeV LHC 
will be presented. 



4 Kinematic distributions for the DDY process 

In this section we study differential cross sections for the production of four-muon final 
states originating from a decay of two 7* after the set of basic cuts has been applied. 
To see the effects of parton showering and intrinsic pt smearing, we present three sets 
of predictions for the SPS and DPS production mechanism: 

• Parton level only (labelled as "PL" in the plots), 

• Parton level + parton shower (labelled "PS"), 

• Parton level + parton shower + intrinsic pt smearing (labelled "PS+cx"). 
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Figure 3: The DPS and SPS single /i distributions in (a) 77 and (b) pt for the process 
pp — )■ 7*7* — at the LHC with 14 TeV. 



The rapidity and transverse momentum distributions for a single // are shown in 
Fig. [3l We see that in the entire 77 and pt range considered here, the SPS background 
dominates over the DPS signal. The r] distribution for the SPS background is more 
central within the allowed rj range. This is because the four muons originate from a 
single hard process, and it is more likely for all of them to be within the 77 acceptance 
if the four-muon system falls in the centre of this allowed region. The pt distributions 
of SPS and DPS muons peak at about the same value, with the SPS distribution being 
slightly harder. Apart from the DPS pt distribution, we see that the other distributions 
shown in Fig. [3] are not sensitive to ISR or intrinsic pt- 

Next we study the separation of DPS and SPS events using their different kinematic 
behaviour in the transverse plane. As discussed in Section [2l at parton level, provided 
that the fi~^fi~ pairs are grouped correctly, the vector pj- sum of the two OS pairs in 
a DPS event will be balanced separately, i.e. Pt = 0. For the SPS, such pt balancing 
is generally achieved only among all 4 muons, which implies that the Pt of the two 
pairs will be opposite to each other, resulting in a back-to-back configuration 
in the transverse plane. These differences can in principle be used to distinguish the 
DPS from the SPS process. In addition, these can be used to find the pairing 
that corresponds to the actual DPS process, and allow further analysis based on the 
kinematic properties of these pairs, i.e. the 7* in our present context. 
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Figure 4: (a) The DPS and SPS differential distributions in the p^-balance variable S 
for the combination that minimises S, and (b) comparison of the DPS differential 

distributions for the "right" , obtained by minimizing S, and "wrong" fi~^fi~ pairs with 
the theoretical distribution obtained by always pairing a fi'^fi~ pair from the same 
7* together. In plot (b), the top (bottom) three keys correspond to parton shower 
simulations without (including) intrinsic Pt smearing. 



To select one of the two combinatoric ways to form two OS muon pairs, we select 
the combination that minimises the following pair-wise balancing variable {S): 

^ _ _ ij^ ^ ij^ \ 0<S <1 4 

2 V PT^^+ + PTt,- Pt4 + PT^,- J 

For the DPS signal, at parton level the correct pairing always leads to 5 = 0, while 
ISR and intrinsic pt effects will distort the picture. Other variables similar to 5* have 
been defined in the literature, for e.g. those in Ref. We have checked that their 
performance are similar to the S variable defined in Eq. H] 

The 5* distributions are displayed in Fig. HI As we can see, the DPS parton level 
result is in strict contrast to the SPS process. In fact, from Fig. it can be seen 
that the values of S close to 0.5-0.7 are the most preferred for SPS, whereas S peaks 
sharply in the lowest bin for the DPS at parton level, as expected. Unfortunately, 
radiation and intrinsic pt effects change the picture dramatically. While the S distri- 
bution for the SPS process gets modified in a relatively moderate way, the additional 
gluon radiation in the DPS process spoils the px balance to a large extent, flattening 
and spreading the distribution up to high values of S. This, of course, significantly 
decreases the power of the variable 5* to discriminate between the SPS and DPS mech- 
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Figure 5: The DPS and SPS differential distributions in the angle A0 between the two 
/i"*" for the process pp — > 7*7* — ?■ fi~^ fi~ fi'^ fi~ at the LHC with 14 TeV. 



anisms. In Fig. Hb we compare the 5* distributions between the pairings obtained in the 
above procedure ("right"), the other pairing ("wrong") and the theoretical distribution 
("theory") obtained by insisting that the fi~^fi~ pair comes from the same 7*. We see 
that the "theory" and "right" distributions are broadly similar, indicating the effec- 
tiveness of the minimum S criterion to pair the OS muons coming from decay of the 
same 7* in DPS, despite the presence of radiation and intrinsic pt effects. We estimate 
that using this requirement leads to selecting the correct fi~^fi~ pair for around 75% 
of the DPS event sample. For the SPS background, due to much more complicated 
kinematical features, the minimum S criterion does not select any particular physical 
SPS configurations. 

We next study the distribution in the angle on the transverse plane between the two 
same-sign muons, A0, here chosen to be /i^/i^, cf. Fig. |5l The SPS events tend to occur 
with bigger angular separation, which can be traced back to the tendency for the 
pairs to be back-to-back in the transverse plane at parton level. For the DPS signal, 
one would expect that independent scatterings lead a fiat A0 distribution. However, 
the requirement that all four /i"*"//" pairs pass the invariant mass cut preferentially 
selects events with all OS /i-pairs apart, resulting in a smaller angular separation for 
same-sign muons. The effect for the opposite-sign muons can be observed in Fig. |H1 
which shows the angular distributions for the "right" and "wrong" combinations of 
pairs. The "right" pairs minimize S, making it more likely for DDY that the two 
paired muons come from decay of the same 7*. In particular, at LO where the DY 7* 
have no pt, choosing S to be minimal tends to pick out only those ^'^ combinations 
which are back-to-back, see Fig. [6^. The preference towards larger angular separation 
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Figure 6: The DPS and SPS differential distributions in the angle A0 between the 
opposite-sign muons for the (a) "right" and (b) "wrong" pairs for the process pp — )■ 
7*7* at the LHC with 14 TeV. 



between "right" DPS pairs remains after including radiation and smearing effects. 
However, the SPS events also tend to have larger angular distributions between the 
OS muons in pairs that minimize the value of 5*. This tendency is no longer true for 
SPS "wrong" fi~^fi~ pairs, i.e. the combination with a larger value of S. In fact, the 
qualitative effect is opposite to that of /i"*"//^ shown in Fig. [5l which is not unexpected 
since both distributions involve azimuthal angular separations between muons from the 
wrong pairs, and would be correlated to certain extent. 

Further differential distributions involving /i^/i" pairs, selected using the S'-variable 
criterion, are shown in Figs. [7] (p^ and rapidity y) and [8](A0 and Ay). As expected, the 
Pt distribution of the DPS muon-pairs is peaked at very low values, lower than those 
for the SPS distribution, even after the showering and smearing effects are included, 
cf. Fig. [7|i. In particular, in the low pt region an excess of events compared to a pure 
SPS model prediction might be observed. As shown in Fig. [7b, both the SPS and DPS 
have fairly symmetric rapidity distributions for the constructed pairs, with very 

similar distributions. In Fig. [8^ we present the distributions for the angle (in the trans- 
verse plane) between the fi~^fi~ pairs from the "right" combination, A(j){fi~^ fi~ , fi'^ fi~) . 
The DPS distribution is relatively uniform over the whole range of A0(/i"'"/i~, 
consistent with the assumption of two independent scatterings in the DDY process. 
The SPS contribution, on the other hand, peaks at A(j){fi~^fj,~ , fj,~^fj,~) = vr, which is a 
consequence of the transverse momentum conservation for the SPS process at leading 
order, as already discussed. Fig. |8]d shows the properties of the SPS and DPS processes 
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Figure 7: The DPS and SPS differential distributions in (a) pt and (b) y for the "right" 
/i"*"/!" pairs for the process pp — )■ 7*7* — )■ fi'^ fi~ fi~ at the LHC with 14 TeV. 
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Figure 8: The DPS and SPS differential distributions in (a) the angle A0 between the 
two pairs and (b) the absolute value of the difference in rapidity \Ay\ for the 

"right" muon-pairs for the process pp — > 7*7* — > fi~^ fi~ fi~^ fi~ at the LHC with 14 TeV. 



with respect to correlations in the longitudinal direction, or more precisely the absolute 
value of the difference between the rapidities of the two pairs, \Ay{fi~^ fi~ , fi~)\. 

As opposed to distributions depending on variables measured in the transverse plane. 
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Figure 9: The DPS and SPS differential distributions for all four muons in (a) the 
invariant mass m(4/x) and (b) the transverse momentum pT(4/i) for the process pp — 
7*7* at the LHC with 14 TeV. 



distributions involving rapidity variables exhibit very little dependence on radiation 
and intrinsic smearing, an effect first noticed and put to use in [27|. We see that as 
\Ay\ increases, the SPS contribution falls more steeply than that of DPS, although the 
difference is not as strong as in the double J/ip analysis in 27 . 

Yet another possibility for studying the DPS and SPS distributions is offered by 
the distributions sensitive to quantities constructed from four vectors of all muons in 
the final state. In Fig. [9] we investigate the invariant mass m(4/i) and the transverse 
momentum ^^(4//) distributions of all four muons. A significant difference in the invari- 
ant mass distribution can be observed for the SPS and DPS mechanisms, cf. Fig. [9^, 
and an excess of DPS events might be observed in the low m(4/i) region. This can 
be traced back to the lower cut on the fj^~^fi~ invariant mass, which, as already noted, 
favours configurations with all OS /x-pairs apart. For DPS the impact of the cut is 
relatively mild, as the independent scattering hypothesis implies that OS pairs from 
different hard scatterings tend to be well separated. This is not the case for the SPS 
background, which results in higher invariant mass of all leptons in the final state. 
From this figure, it can be inferred that events with m(4/i) > 40 GeV are dominated 
by the SPS background, and hence could form a control region for the SPS background 
estimation. A more aggressive cut on max m(4/i), down to around 30 GeV, appears 
to be possible. Such a cut would reduce the SPS background by around 30%, while 
having a relatively mild impact on the signal. However, we remain conservative here 
and refrain from doing so. The pT(4/i) distribution shown in Fig. [9]d measures the 
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(a) DPS and (b) SPS mechanism for the process pp 7*7* — j- jX^ jX^ jjT at the LHC 
with 14 TeV. 



sensitivity of the SPS and DPS to extra radiation, simulated by parton showering, 
and intrinsic smearing. Both distributions are affected by these effects in a similar 
way, with the peak of the showered distribution at approximately the same value and 
moving to higher px as the smearing is switched on. 

Since the two types of correlations, i.e. in the transverse plane and in the longitudi- 
nal directions, between the muons in the final state probe different kind of information 
on the same process, it is interesting to study the double differential distributions in 
A0(/i"'"/x~, /i"*"/!") and |A?/(/i"'"yU~, yU"*"/!")!. In Fig. [10] we present scatter plots in the 
A0(yU~'"/i~, /i^/i") and |A?/(/i+/i~, yU+zi")! plane for the DPS and SPS samples. The 
former appears to be slightly more populated in the region with small A0 and |A|/|, 
whereas events of the latter tend to occupy the corner of large A0(/i^/i~, and 
small I A|/(/i"'"/i~, /i"*"//") |. 

Concluding, observing DPS through four-muon final states from DDY scatterings 
seems to be a challenging task due to its low production rate, cf. Table [H Although an 
excess of DPS events could be observed in the low pt{^^ and m(4/i) regions, the 
studied differential distributions between the SPS and DPS event samples do not show 
kinematical differences that are striking enough to justify putting additional hard cuts 
(perhaps with the exception of m(4/i)) at the cost of further reducing the signal yield. 

On the other hand, we see that in addition to pt{^X^ l^iT) and m{4fi), variables such 
as S, A(f){fi~^fi~,fi^fi~), |A?/(/i"'"/i~, /i^/i")! or A0(/i+/i~) of the wrong muon-pairs also 
show moderate, non-trivial differences between the DPS and SPS components, and we 
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are able to understand the origin of their different distributions. To fully utilize these 
differences, it would be more effective to apply the template method, as is done for 
example in the 7 + 3j studies in CDF and DO to extract the DPS fraction. Moreover, 
we see from the double differential distribution in Fig. (TU] that additional information 
can be obtained by simultaneously considering more than one kinematic quantity. 



5 Kinematic distributions for double J/ijj produc- 
tion 

As seen in Section [3l after applying sets of basic cuts the total cross sections for double 
J / ip production at LHCb are significantly higher than the equivalent cross sections for 
the DDY process. Even more importantly the S/B ratio is also better. This makes 



double J/ijj production an excellent candidate process to observe DPS |27H29 . 

Methods to study double J/ip and to extract the DPS signal from the SPS back- 
ground have been discussed in some detail in Here we extend this study by 
comparing the different properties between double J/ip and DDY in pp collisions at 14 
TeV. 

There are a few important differences between double J/ip and DDY. First, their 
production is initiated by different partons: while double J/ip production is driven by 
the gluon-gluon fusion process, both SPS DDY and the single DY subprocess in DPS 
DDY require, at lowest order, a quark and an antiquark. The longitudinal momentum 
densities of these partons are different, which would lead to different rapidity distribu- 
tions for the pairs. Different initial-state partons also affect the parton shower, 
and hence pt- However, the hard scales of the double J/ip and DDY are also different, 
which will in turn affect the amount of QCD radiation. The combined effect has an 
impact on pt as well as azimuthal angular separations in the transverse plane of the 
reconstructed objects^ 

In Fig.[TT]we show the rapidity and A0(/i"'"/i~, /U.+/i~) distributions of the /i^/i" pairs 
from double J/ip production. Here the invariant mass of a pair readily determines 
if they originate from a J/ip resonance or not, and hence there is no ambiguity in pairing 
up the OS muons. Compared with DDY (c/. Fig. Ub), we see in Fig. [TTk that both the 
double J/ip SPS and DPS are more central, with the effect of SPS stronger. This is 
due to the smaller gluon densities at high parton x compared to the light (anti) quarks, 
making high 77 values for the iJ,'^fi~ pairs less likely. In Fig. [TTb the impact of QCD 
radiation on A0(yU"''yU~, /x^/x") is clearly seen, where the distribution including parton 
shower peaks slightly towards zero, instead of vr, the parton level result. This should 
be compared with the corresponding figure for DDY, Fig. [8^, which shows a milder 
radiation effect. 



^In the present study we do not attempt to disentangle the effects on pT from these sources. 
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Figure 11: The DPS and SPS differential distributions in (a) y and (b) A0(/i+/i , /i+yU ) 
of pairs for the process pp J /ipJ /ip — )■ jX^ iJi~ jX^ n~ at the LHC with 14 TeV. 



A consequence of the well defined way to pair up the OS muons is that it is not 
necessary to introduce the variable 5* as a way to group the muons, although it could 
still be used to help distinguish the SPS and DPS components. Since for the SPS DDY 
process the pairing according to value of S does not correspond to finding two /I's from 
a decay of the same 7*, more significant differences between DDY and double J /ip are 
expected for the SPS background. 

The four-muon invariant mass distribution for double J/ip production is shown in 
Fig.[T2b. The DPS component peaks at higher value of m(4/x) than the SPS component. 
The reverse is observed for DDY, c/. Fig. [9^. This can partly be attributed to the /i"*"//" 
invariant mass cut applied in the DDY analysis, which will tend to push the SPS m(4/i) 
distribution towards higher values. Of course there is no equivalent mass cut in the 
J / study, since there is no ambiguity in pairing up the muons. A quantity related to 
m(4//) is |A|/(/i+yU~, /i+yU")!, as in the centre of mass frame, a larger m(4/i) will implies 
a larger |A|/(/x"'"/i~, given a fixed pt(/^^/^~)- This distribution is displayed in 
Fig. [T2b . The observation for double J/ip oi a more steeply falling \Ay\ for the SPS 
is consistent with the more steeply falling m(4/i). For DDY, the broader SPS m(4/x) 
distribution implies that \Ay\ would be more slowly falling, thus making \Ay\ less 
effective as a variable to extract DPS from SPS component in DDY. 

Finally, we note that overall our current studies reinforce the statements made 



m 



27j regarding the importance of radiation and intrinsic pt smearing effects for the 



observables involving quantities measured in the transverse plane. At the same time, we 
confirm the relative insensitivity to these effects on the variables measuring correlations 
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Figure 12: The DPS and SPS differential distributions for (a) the invariant mass m(4/i) 
and (b) the absolute value of the difference in rapidity |A|/(/i"'"/x~, /i^/i")! for the fj.~^f^~ 
pairs for the process pp — )■ J/ipJ/ili — )■ ii^ ii~ n'^ ii~ at the LHC with 14 TeV. 



in the longitudinal direction. Therefore it could be beneficial to investigate in more 
detail how "longitudinal variables" may be used to facilitate other DPS searches in the 
future. 



6 Conclusions and outlook 

Measurements undertaken at the LHC will be crucial for improving the theoretical 
description of multiple parton scattering. In this paper we have studied the production 
of two OS muon pairs via two hard scatterings of the DY type. DY production is 
considered a standard candle process at hadron colliders. The leptonic final states 
provide very clean experimental signatures, while the status of theoretical predictions is 
under control. All this makes the production of two OS muon-pairs a model candidate 
process for probing the nature of DPS. 

Depending on the invariant mass of the OS muon pairs, they may originate either 
from continuum (7*) or a hadronic resonance, e.g. J/ip. In the course of our study, 
we have found that at 7 TeV, the DDY cross section and luminosity are too low for 
this process to be observed with the LHCb detector. The prospects for observing 
DDY should improve at 14 TeV LHC given the increases in cross section, S/B ratio 
and luminosity. On the other hand, the production of the same final state through a 
double J /ip process possesses characteristics which should allow a measurement of DPS 
at the LHCb in pp collisions at both 7 and 14 TeV. However, we stress that these two 



19 



processes are sensitive to different initial state partons ((anti)quark-(anti)quarlc for 
DDY and gluon-gluon for double J/ip)- Hence they probe different correlation effects 
and provide complementary input to double parton distributions. 

For the same reason, it might prove valuable to additionally investigate four muon 
final states where one muon pair comes from a decay of 7* while the second originates 
from the decay of a low-mass resonance, e.g. J The cross section for such a process 
might be expected to lie between double J/ip and DDY, and hence is likely to be 
observable. Furthermore, although in this paper we have concentrated on pairs of 
J/^/^'s as resonant sources of four-muon final states at low invariant masses, another 
promising opportunity to study DPS can be offered by studying higher mass resonances, 
in particular involving T production. For example, measurement of a^s for both double 
J / ip and double T processes could provide information on the scale dependence of the 
two-gluon parton distribution function. This would complement well the study of 
scale dependence that could be performed in e.g. 7 + 3j, H^^H^^ and other jet-based 



signatures proposed for the general purpose detectors ATLAS and CMS [39 

Up on completion of this paper, a suggestion to study double T has been put forward 

in 



r 
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